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Abstract 

We study the Poissonnian ensembles of Markov loops and the associated renor- 
malized self intersection local times. 



1 Introduction 



The purpose of this paper is to explore some simple relations between Markovian path and 
loop measures, spanning trees, determinants, and Markov fields such as the free field. The 
main emphasis is put on the study of occupation fields defined by Poissonian ensembles 
of Markov loops. These were defined in [9] for planar Brownian motion in relation with 
SLE processes and in [10] for simple random walks. They appeared informally already 
in [24]. For half integral values | of the intensity parameter a, these occupation fields 
can be identified with the sum of squares of k copies of the associated free field (i.e. 
the Gaussian field whose covariance is given by the Green function). This is related to 
Dynkin's isomorphism (cf [2], [17], [13]). We first present the results in the elementary 
framework of symmetric Markov chains on a finite space, proving also en passant several 
interesting results such as the relation between loop ensembles and spanning trees. Then 
we show some results can be extended to more general Markov processes. There are no 
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essential difficulties when points are not polar but other cases are more problematic. As 
for the square of the free field, cases for which the Green function is Hilbert Schmidt such 
as two and three dimensional Brownian motion can be dealt with through appropriate 
renormalizat ion . 

We can show that the renormalised powers of the occupation field (i.e. the self inter- 
section local times of the loop ensemble) converge in the two dimensional case and that 
they can be identified with higher even Wick powers of the free field when a is a half 
integer. 

2 Symmetric Markov processes on finite spaces 

Notations: Functions and measures on finite (or countable) spaces are often denoted as 
vectors and covectors. 

The multiplication operator defined by a function / acting on functions or on measures 
is in general simply denoted by /. but sometimes it will be denoted Mf. The function 
obtained as the density of a measure /i with respect to some other measure u is simply 
denoted -. 

Our basic object will be a finite space X and a set of non negative conductances 
Cx,y — Cy^x, indexed by pairs of distinct points of X. 

We say {x, y} is a link or an edge iff C^^y > and an oriented edge {x, y) is defined by 
the choice of an ordering in an edge. We set —{x, y) — {y, x) and if e = {x, y), we denote 
it also (e~, e"*"). 

The points of X together with the set of non oriented edges E define a graph. (X, E). 
We assume it is connected. The set of oriented edges is denoted E". 

An important example is the case in which conductances are equal to zero or one. 
Then the conductance matrix is the adjacency matrix of the graph: Cx,y = l{x,y}£E 

2.1 Energy 

Let us consider a nonnegative function k on X. Set \x = f^x + Xlj/ ^x,y Py = P is 

a A-symmetric (sub) stochastic transition matrix: )^xPy = \Px with = for all x in 
X and it defines a symmetric irreducible Markov chain 

We can define above it a A-symmetric irreducible Markov chain in continuous time Xt, 
with exponential holding times,of parameter 1. We have Xt — ^Nt-i where Nt denotes a 
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Poisson process of intensity l.The infinitesimal generator writes Ly — Py — 5y. 

We denote by Pt its (sub) Markovian semigroup exp{Lt) — L and Pt are 

A-symmetric. 

Wc will consider the Markov chain associated with C, k, sometimes in discrete time, 
sometimes in continuous time (with exponential holding times). 
Recall that for any complex function z^,x & X, the "energy" 

e{z) = {-Lz,z)^ = ^ -{LzYz'^K 

is nonnegative as it can be written 

e{z) = ^ J] C,,y{z- - zy){z^ - ^) + J] i^xz""^ = - E 

x,y X X x,y 

The Dirichlet space ([4]) is the space of real functions equipped with the energy scalar 
product defined by polarization of e. 

Note that the non negative symmetric "conductance matrix" C and the non negative 
equilibrium or "killing" (or "equilibrium") measure k are the free parameters of the model. 

We have a dichotomy between: 

- the recurrent case where is the lowest eigenvalue of —L, and the corresponding 
eigenspace is formed by constants. Equivalently, PI — 1 and k vanishes. 

- the transient case where the lowest eigenvalue is positive which means there is a 
"Poincare inequahty": For some positive £, the energy e(/, /) dominates £ (/, /);^ 
for all /. Equivalently, k does not vanish. 

Wc will now work in the transient case. Wc denote by V the associated potential 
operator (— L)^^ = Pfdt. It can be expressed in terms of the spectral resolution of L. 

We denote by G the Green function defined on as G"*'^ = y" = ~ P)~^]y 

i.e. G = (Ma — G)~^. It induces a linear bijection from measures into functions. We set 

(G/.)^ = E,G-'>. 

Note that e{f,GiJi) = {f,ix) (i.e. J^xf'^f^^) function / and measure In 

particular = 1 as e(l, /) = ^ fi^x = (/, !)«• 

See ([4]) for a development of this theory in a more general setting. 



3 



In the recurrent case, the potential operator V operates on the space A of functions / 
such that (/, = as the inverse of the restriction of J — P to A-'-. The Green operator G 
maps the space of measures of total charge zero onto A"*-. Setting for any signed measure 
u of total charge zero Gu = Vj. we have for any function /, {v, f) = e{Gv, f) (as 
e{Gy, 1) = 0) and in particular — P = e{G{Sx — Sy), /). 

2.2 Feynman-Kac formula 

For the continuous time Markov chain xt (with exponential holding times) and k{x) any 
non negative function, we have the Feynman Kac formula: 



For any nonnegative measure X; set = (— L + ) ^ and G^ = V^Mi = {Mx + 
— C)^^. It is a symmetric nonnegative function on X x X. Go is the Green function 

1 1 2 

G, and G^ can be viewed as the Green function of the energy form e^^ = e + || |Il2(x)- 

Note that has the same conductances C as e, but x is added to the killing measure. 
Note also that V-^ is not the potential of the Markov chain associated with when one 
takes exponential holding times of parameter 1 but the Green function is intrinsic i.e. 
invariant under a change of time scale. Still, we have by Feynman Kac formula 



2.3 Countable spaces 

The assumption of finitcncss of X can be relaxed. On countable spaces, the previous 
results extend easily when under spectral gap conditions. In the transient case, the 
Dirichlet space H is the space of all functions / with finite energy e(/) which are limits 
in energy norm of functions with finite support. The energy of a measure is defined as 
sup/gH^^^- It includes Dirac measures. The potential G/j, is well defined in H for all 
finite energy measures /i, by the identity e(/, G/i) — (/, /i), valid for all / in the Dirichlet 
space. 

Most important cases are the non ramified covering of finite graphs. 



= [exp(i(L - M,)] 




We have also the "resolvent" equation V -V^^^ VMkV^ = V^MkV . Then, 

G-G^^ GM^G^ = G^M^G 
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3 Loop measures 

3.1 A measure on based loops 

We denote the family of probability laws on piecewise constant paths defined by Pt- 

IP"'(7(^i) = xi, ...,7(tfe) =Xh) = Pt^{x,xi)Pt^-t^{xi,X2) . . .Pt^-t^_-,ixh-i,Xh) 
Denoting by ^(7) the number of jumps and Tj the jump times, we have: 

^x{p{l) = k, 7ti = xi, 7t;,_i = Xk-i, Ti e dti, E dtk) 



^X^X2---^Xk 

For any integer p > 2, let us define a based loop with p points in X as a couple I — r) = 
((Cm, 1 < "T. < p), (Tjn, 1 < m < p + 1), ) in XP X R^"^^, and set — ^p+i (equivalently, 
we can parametrize the the discrete based loop by Z/pZ). The integer p represents the 
number of points in the discrete based loop ^ = (Cij will be denoted p(0- 

Note two time parameters are attached to the base point since the based loops do not in 
general end or start with a jump. 

Based loops with one point (p = 1) are simply given by a pair r) in X x ]R_,_. 

Based loops have a natural time parametrization l{t) and a time period T{^) = 

If we denote YnLin by T^: /(t) = Cm-i on [T^_i,T^) (with by convention 

ro = Oandeo = y- 

A (7-finite measure /i is defined on based loops by 



/CO ^ 



where P^'^ denotes the (non normalized) "law" of a path from x to y oi duration t : If 

ti < t2... <th<t, 

¥^'y{l{t,)^x^,..J{t,)^x,) = [PtXAPt2-tXl-4Pt-t,r/Y 

Ay 

Its mass is p^'^ = . And for any measuable set A of piecewise constant paths indexed 
by [0 t] , we can also write 

r,^'{A)=Fx{An{xt = y})^ 

Ay 



Prom the first expression, we see that by definition of if <t2--- < th <t, 

= x^, ^xn,Te dt) = [P*,+,_,J^jP,,_*J^^..[P*,_*,_J^-i^cit (1) 

Note also that for /c > 1, using the second expression of P^'^ and the fact that condition- 
ally to Nt = k, the jump times are distributed like an increasingly reordered A;— uniform 
sample of [0 t] 

A^P,^'^(p = k, ^2 = X2. ^ Xk,Ti e dh, ...,Tk e dtk) 
^P^^P^^...P^n^o<t,<...t,<t}e-'dh...dh 

Therefore 

Kp = ^, 6 = xi, Ck = Xk, Ti e dh, .., Tk e dtk, T e dt) (2) 
^ P;^.P;f ^^'"'''^■■■''''''^'K -'dh...dtkdt (3) 

for k > 1. 

Moreover, for one point-loops, — — ^i: "^i £ dt} — ^dt 

3.2 First properties 

Note that the loop measure is invariant under time reversal. 

If D is a subset of X, the restriction of to loops contained in D, denoted /j,^ is 
clearly the loop measure induced by the Markov chain killed at the exit of D. This can 
be called the restriction property. 

Let us recall that this killed Markov chain is defined by the restriction of A to D and 
the restriction of P to (or equivalently by the restriction cd of the Dirichlet norm 
e to functions vanishing outside D). 

As J ^-^e~*dt = |, it follows from (2) that for A; > 1, on based loops, 

Kp(0 = ^,6 = ^1, ...,^k = X,) = Ip^^.-.P^I: (4) 
In particular, we obtain that, for A; > 2 

^{p = k) = ]-Tr{P'') 
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and therefore, as Tr{P) — 0, 



liip >l)^Y. rTriP') = - log(det(7 - P)) = log(det(G) J] A.) 

2 X 

since (denoting Mx the diagonal matrix with entries A-^), we have 

^ ' det(Mj) 

Moreover 

/oo 
p(/)l{,>i}/.(d/) = J^TrlP'^) = Tr((/ - P)-ip) = Tr{GC) 
2 

3.3 Loops and pointed loops 

It is clear on formula 1 that /i is invariant under the time shift that acts naturally on 
based loops. 

A loop is defined as an equivalence class of based loops for this shift. Therefore, n 
induces a measure on loops also denoted by /i. 

A loop is defined by the discrete loop ^° formed by the in circular order, (i.e. up to 
translation) and the associated scaled holding times. We clearly have: 

//(e° = (Xi,X2,...,X,)°) = P,^^..P,^^'= 

However, loops are not easy to parametrize, that is why we will work mostly with 
based loops or pointed loops. These are defined as based loops ending with a jump, or as 
loops with a starting point. They can be parametrized by a based discrete loop and by 
the holding times at each point. Calculations are easier if we work with based or pointed 
loops, even though we will deal only with functions independent of the base point. 

The parameters of the pointed loop naturally associated with a based loop are Ci, ...,Cp 
and 

Ti + Tp+i= T*, Ti^T*, 2<i<p 

An elementary change of variables, shows the expression of /i on pointed loops writes: 
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Ii(p ^k,^i^ Xi,r* e dt,) = P^l...P^^^e-^'^dh...dtk (5) 

Trivial (p = 1) pointed loops and trivial based loops coincide. 
Note that loop functionals can be written 

with $fc invariant under circular permutation of the variables (Ci,T*). 
Then, for non negative 

and by invariance under circular permutation, the term ti can be replaced by any ti. 
Therefore, adding up and dividing by A;, we get that 

j ^k{l°Mdl)^ J^J2^k{xiA)P^^-P^,'e-^''dt,...dtk 

The expression on the right side, applied to any pointed loop functional defines a 
different measure on pointed loops, we will denote by fi* . It induces the same measure as 
11 on loops. 

We see on this expression that conditionally to the discrete loop, the holding times of 
the loop are independent exponential variables. 

^Jnp = K^i = x,,T* edt,)^^ n c^^,^^^^e-'^du (6) 

Conditionally to p(^) = A;, T is a gamma variable of density -j^^e"* on IR+ and 

1 < i < k) an independent ordered fc-sample of the uniform distribution on (0,T) 
(whence the factor j). Both are independent, conditionally to p of the discrete loop. We 
see that /x, on based loops, is obtained from ^ on the loops by choosing the based point 
uniformly. On the other hand, it induces a choice of biased by the size of the r*'s, 
different of //* (whence the factor |. But we will consider only loop functionals. 
It will be convenient to rescale the holding time at each by A^. and set % — 
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The discrete part of the loop is the most important, though we will see that to estab- 
lish a connection with Gaussian fields it is necessary to consider occupation times. The 
simplest variables are the number of jumps from x to y, defined for every oriented edge 

Nx,y = #{i ■^i = x, = y} 
(recall the convention ^p+i = ^i) and 

y 

Note that — #{i >^'-ii = x} except for trivial one point loops for which it vanishes. 
Then, the measure on pointed loops (5) can be rewritten as: 

yri*(p = = a;,r e (it) = e"^"*y and (7) 



x,y X ieZ/pl 

Another bridge measure //^'^ can be defined on paths 7 from x to y: ^^''^{d'y) — 

Note that the mass of ^^'^ is G^'^ . We also have, with similar notations as the one defined 
for loops, p denoting the number of jumps 

^f^y{p{-f) = k, = xi, 7t,_i = Ti e dh, Tk-i G dtk-i, T e dt) 

Cx,X2Cx2,X3---Cx^._i,y „-tj+ J+ JJ. 

— r—: r ^{o<ti<...<tk<t}^ aii...atkai 

'^x'^xi-'-'^y 

3.4 Occupation field 

To each loop 1° we associate local times, i.e. an occupation field {lx-,x G X} defined by 

fT(0 -, p(0 



Jo h{s) 



for any representative I = {^i, r*) of 1°. 

For a path 7, 7 is defined in the same way. 
Note that 



M(l-e-"^^)l{,=i})= re-*(l-e-^*)^ = log(l + f ) (9) 

Jo -^x 
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(by expanding 1 — e~^* before the integration, assuming first a small and then by ana- 
lyticity of both members, or more elegantly, noticing that J^ie"'^ — e"*^^)^ is symmetric 
in (a, b) and (c, d)). 
In particular, n{l^l^p=ij) = j-. 

From formula 5, we get easily that the joint conditional distribution of {l^, x G X) 
given {Nx, x G X) is a product of gamma distributions. In particular, from the expression 
of the moments of a gamma distribution, wee get that for any function $ of the discrete 
loop and A; > 1, 



K(iyMp>i}^) = K + k- 1)...(N, + 1)7V,$) 

In particular, ^(P) = -^[fiiN^) + 1] = G'"'''. 

Note that functions of / are not the only functions naturally defined on the loops. 
Other such variables of interest are, forn > 2, the multiple local times, defined as follows: 

- 1 

lXU...,Xr. = ^ / l{^{tl)=Xl+J,....at„-j)=Xr.,...atu)=Xj} H ^^^^ 

j=0 J0<ti<...<t„<T ^Xi 

It is easy to check that, when the points Xi are distinct, 

n—l n 

^"'•■■"" = E E nite.=W^r (10) 

j=0 l<ii<..<in<p(l) 1=1 

Note that in general pi'-.^^fc cannot be expressed in terms of /. 

li xi = X2 — . . . — Xn, pi'-'^" = (n-i)! [^^]"- vicwcd as a n-th self intersection 

local time. 

One can deduce from the defintions of the following: 

Proposition 1 ;^(pl.-,=^n) ^ QXuX2QX2,X3 __ _QXn,X^ 

Proof. Let us denote j- [Pt]y by p^'^ or pt{x,y). From the definition of pi' - .^n and 
fjL, - equals: 

X]A^^ / / ^pti{x,xi+j) . . .pt-t„{xn+j,x)Yldtidt 

X j=0 J{0<ti...<t„<t} ^ 
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where sums of indices k + j are computed mod(n). By the semigroup property, it equals 



n-l 



^ TPt2-ti {Xl+j, X2+j) . . . Pti+t-t„ {Xn+j, Xl+j) Yl dtidt. 

j=0 J{0<ti<...<t„<t} I 

Performing the change of variables V2 — t2 — ti, ..,Vn — tn — tn-i, Vi — ti + t — tn, and 
V — ti, we obtain: 

n-l „ ^ 

52 / ] ] Pv2 {Xl+j, X2+j) ...Pvi {Xn+j, Xl+j) TT dVidv 

~^J{0<v<vi,0<Vi}Vl + ---+Vn 

n—l p 

= / . V Pv2{xi+j,X2+j) Pv,{Xn+j,Xi+j)T\dVi 

j^O J{0<Vi} Vi + ...+Vn 

n „ 

= y^l Py^{Xi,X2) . . .Pv^{Xn,Xi)Y\dVi 

j^^J{0<Vi}Vl + ...+Vn 

J{0<Vi} 

_ QXl,X2QX2,X3 QXn,Xl 

Note that another proof can be derived from formula (10) ■ 

Let us come back to the occupation field to compute its Laplace transform. Prom 
the Peynman-Kac formula, it comes easily that, denoting Mk the diagonal matrix with 
coefficients ^ 

Pr(e-(^>^) - 1) = f (exp(i(P -I-Mjfx- exp(i(P - /))^). 
Integrating in t after expanding, we get from the definition of /j, (first for x small enough): 

(e-(''>^) - l)d/.(/) = Y,J^ [Tr{{P - M^f) - Tr{{Pf)]^e-'dt 



00 

Y.T[Tr{{P-M^)')-Tr{{Py 



k 

k=l 



-Tr(log(/-P + M^)) +Tr(log(/-P)) 
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Hence, as T'r(log) = log(det) 

J (e-('^^^) - l)dii{l) - log[det(-L(-L + M^/a)"^)] = -logdet(/ + VM, 



which now holds for all non negative x both members are analytic in x- Besides, by 
the "resolvent" equation: 



det(7 + GM,)-^ = det(7 - G^M^) = (11) 

Note that det(/+G'M^) = det (J+M^GM^^) and det(/-G^M;^) = det(/-M^G;^M^), 
so we can deal with symmetric matrices. Finally we have the 

Proposition 2 /.(e-i^'^) - 1) = - log(det(7 + M^GM^)) = log(^) 

Note that in particular //(e"*'"" — 1) = — log(l + tG^'^). 
Note finally that if x bas support in D, by the restriction property 

/^(l{W)=o}(e"^''^ - 1)) = -log(det(/ + M^G^M^)) = log(^^^) 

Here the determinants are taken on matrices indexed by D and G^ the Green function 
of the process killed on leaving D. 

For paths we have P^'^(e"(''^)) = ^ exp(t(L - M^))^^y. Hence 

Ay 

Also E^(e-<^''^>) = ^ylG^j'^'^Ky i.e.[G^K]''. 

4 Poisson process of loops 
4.1 Definition 

Still following the idea of [9], which was already implicitly in germ in [24], define, for all 
positive a, the Poissonian ensemble of loops Ca with intensity a/i. We denote by P or 
its distribution. 
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Recall it means that for any functional $ on the loop space, vanishing on loops of 
arbitrary small length, 



^(e^E(e£a*(0 = exp(a J (e**« - 



Note that by the restriction property, = {/ G Ca,l C D} is a Poisson process of 
loops with intensity /i^, and that is independent of Ca\C^- 

We denote by VC^ the set of non trivial discrete loops in Ca- Then, 



with A^i"^-* = -^a;(0 ^^<i ^x^y — Ylfi&Ca ^x,y{l), when these loops are distinct. 

We can associate to Ca a a-finite measure (in fact as we will see, finite when X is 
finite, and more generally if G is trace class) called local time or occupation field 

leCc. 

Then, for any non-negative measure x on X 

E(e-(^'^)) = exp(Q; J (e-('^^) - l)di^{l)) 
and therefore by proposition 2 we have 

Corollary 3 E(e-(^«'>^)) = det(7 + M^GM^)-" = (^g^)" 
Many calculations follow from this result. 

It follows that E(e^*^" ) = (l+tG^'^)^". We see that Ca follows a gamma distribution 
r{a,G^'^), with density l{a;>o}^r^^ (in particular, an exponential distribution of 

mean G^'^ for a = 1). When we let a vary as a time parameter, we get a family of gamma 
subordinators, which can be called a "multivariate gamma subordinator" . 

We check in particular that E(£q ) — aG^'^ which follows directly from /x(/x) — G^'^. 
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Note also that for a > 1, 
More generally, for two points: 



X y 

This allows to compute the joint density of Ca and Ca in terms of Bessel and Struve 
functions. 

We can condition the loops by the set of associated non trivial discrete loop by using 
the restricted cr-field a{VCa) which contains the variables Nx,y. We see from 9 and 7 that 



~r ~K.x 



The distribution of {N^^ , x £ X} follows easily, from corollary 3 in terms of generating 
functions: 

nXiisf'*') = det(5,, + / AA(1-''.)(1-V g,„)- (12) 

X V ^xOy 

SO that the vector of components follows a multivariate negative binomial distribution 

(see for example [26]). 

It follows in particular that A^i""* follows a negative binomial distribution of parameters 
—a and . L^ . Note that for a = 1, Nx"* + 1 follows a geometric distribution of parameter 



KG 
1 



4.2 Moments and polynomials of the occupation field 



It is easy to check (and well known from the properties of the gamma distributions) that 

the factorial moments of A^"^ : 

(iyi"UA-)(A;^V A- -!)...(# + 1) 



the moments of La are related to the factorial moments of A^"^ : 
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It is well known that Laguerre polynomials L^^ with generating function 

oo 



ut 

e 1-* 



ii-ty 





,,a. — l^ — u 



are orthogonal for the T{a,l) distribution with density " p^J^ — l{u>o}- They have mean 
zero and variance Hence if we set = G=^'=^and P^'^ix) = {-(TfLf-^\l), the 

random variables P^''^''{Ca ) are orthogonal with mean and variance (j^'^^^^^j^, for 
k>0. _ _ _ 

Note that Pf''^''(£a ) = Ca ~ (^^x = C.a — IE('Cq. ). It will be denoted ■ 

. . ut 

Moreover, we have T^o^'Pn^) = E(-^^)'=4""'^(^) = ^ 
Note that 

g l + CT^t g 1 + 

^^(l + (7^i)°(l + (7j,s)"'' 

^ :((! - ^^)(1 - i ^((c;-.2/)2)-a 



= (1 - 5^(6-^'^) 2)-". 

Therefore, we get, by developping in entire series in (s, t) and identifying the coefficients: 

E(pr^(£;^), prn>c;')) = 4,(G'-.^)2^^i^±iH^±i:iii (13) 

Let us stress the fact that G^'^ and G^'^ do not appear on the right side of this 
this formula. This is quite important from the renormalisation point of view, as we will 
consider in the last section the two dimensional Brownian motion for which the Green 
function diverges on the diagonal. 

More generally one can prove similar formulas for products of higher order. 
Note that since G^M^ is a contraction, from determinant expansions given in [25] and 
], we have 



det(7 + M^GM^)-" = 1 + X^^...X^,Pera{G.„^^, l<l,m<k) 

k=l 

and then, from corollary 3, it comes that: 

^{(^a,x) ) = Yxh--XikP(^ra{Gi^,i^,l <l,m<k) 
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Here the a-permanent Per^ is defined as X^o-eSfc ■■■^ifc,V(fc) with m{a) denoting 

the number of cycles in a. 

Note that from this determinant expansion follows directly (see [26]) an explicit form 
for the multivariate negative binomial distribution, and therefore, a series expansion for 
the density of the multivariate gamma distribution. 

It is actually not difficult to give a direct proof of this result. Thus, the Poisson process 
of loops provides a natural probabilistic proof and interpretation of this combinatorial 
identity (see [26] for an historical view of the subject). 

We can show in fact that: 

Proposition 4 For any (ii, ...ik) in X", E{Ca ■■■J^a ) = Pera{G''^''"' , 1 <l,m < k) 

Proof. The cycles of the permutations in the expression of PeVa are associated 
with point configTuations on loops. We obtain the result by summing the contributions 
of all possible partitions of the points ii...ik into a finite set of distinct loops. We can 
then decompose again the expression according to ordering of points on each loop. We 
can conclude by using the formula ^{l^^'-'^™-) = G^'^'^'^G^^'^'^ ...G^'^'^'^ and the following 
property of Poisson measures (Cf formula 3-13 in [6]): For any system of non negative 
loop functionals Fj 

■ 

Remark 5 We can actually check this formula in the special case ii = 12 = ■■■ = ik = x. 
From the moments of the Gamma distribution, we have that E((>Cq, )"■) = (G^'^)"q;(q; + 
l)...(a + n — 1) and the a-permanent writes d{n, k)a^ where the coefficients d{n, k) are 
the numbers of n— permutations with k cycles (Stirling numbers of the first kind). One 
checks that d{n + 1, /c) = nd{n, k) + d(n, k — 1). 

Let Si be the set of permutations of k elements without fixed point. They correspond 
to configurations without isolated point. 

Set Per°(G^''^'", l<l,m<k)^ Y.aes° a^^^^G'^^'-w ...G'^Mk) . Then an easy calcula- 
tion shows that: 
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Corollary 6 m{C^\..C^'') = PerUG'^^'^, l<l,m<k) 
Proof. Indeed, the expectation writes 

J2 J2 (-l)*""" n G'^'^'PeraiG'"''', a,bel) 

p<k 7C{l,...ifc},|/|=p lei" 

and 

Pera(G'-''\ a,beI)^J2 H G^''^^er°(G'^"■^^ a,b E J). 

JQi jei\J 

Then, expressing E(£q, ..Xq in terms of Perm's, we see that if J C {1, ...k}, \ J\ < k, the 
coefficient of Per°(G^-^^ a,b e J) is (-l)*'"''^' H ^^^^^ vanishes as (-l)-l-^l = 

(-1)1^1 = and E/Dj(-1)''^^' = (1 - 1)^"'^' = 0. ■ 

Set Qf'iu) = P^'^iu + acr) so that P^'^iC^) = Q^'^i^J)- This quantity will be 
called the n-th renormalized self intersection local time or the n-th renormalized power 
of the occupation field and denoted 

From the recurrence relation of Laguerre polynomials 

we get that 

In particular Q2''^{u) = — 2au — aa'^). 
We have also, from (13) 

E(Qr (Z;^), Qr^iLj)) = + + (14) 

The comparison of the identity (14) and corollary 6 yields a combinatorial result 
which will be fundamental in the renormalizing procedure presented in the last section. 

The identity (14) can be considered as a polynomial identity in the variables a^, Cy 
and C^'K 

If (5fe''^'^(«) = Ylm=o Qm^'^^^x~^-> if denote Nn,m,r,p the number of ordered configu- 
rations of n black points and m red points on r non trivial oriented cycles, such that only 
2p links are between red and black points, we have 
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r p<inf(m,n) 

and therefore 

E E E «'C'?n''^n,m,r-,p = unless p = Z = A;. (15) 

r p<m<k p<n<l 

aQk Qk Nk,k,r,k = (16) 

r 

Note that one can check directly that q'^''^ — ^, and Nk,k,i,k — k\{k — Nk,k,k,k — k\ 
which confirms the identity (16) above. 

4.3 Hitting probabilities 

Let [H^]y = ¥x{xtp = y) be the hitting distribution of F by the Markov chain starting 
at X. Set D ^ F'' and denote , = P)\dxd, = [(/ - P^)]-^ and = 
[{Mx — C)\dxd]~^ the energy, the transtion matrix, the potential and the Green function 
of the process killed at the hitting of F. Recall that 

[H^Ty = + 17 i:.eD[{P''f]lPy = + T..eD\V'']lP'y. Moreover we 

have by the strong Markov property, V = + H^V and therefore G = G^ + H^G. 
(Here we extend and G^ to X x X by adding zero entries outside D x D). 

As G and G^ are symmetric, we have [H^G]y — [H^G]^ so that for any measure u, 
H^{Giy) = G{yH^). 

Therefore we see that for any function / and measure v, e{H^f, G^v) = e{H^f, Gv) — 
e{H^f, H^Gv) = {H^'f, v) - eiH^'f, GiH^'u)) = as {H^'f = 
Equivalently, we have the following: 

Proposition 7 For any g vanishing on F, e{H^f,g) = so that I — is the e- 
orthogonal projection on the space of functions supported in D. 

For further developments see for example ( [12]) and its references. 
The restriction property holds for Ca as it holds for ji. The set of loops inside D is 
associated with /x^ and independent of — C^- Therefore, we see from corollary 3 that 

-(^„-£?,x>^ . det(G^,)det(G^) 
^^'^ ^ = Uet(G) det(G^)^ " 
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Prom the support of the of the Gamma distribution, we see that n(l{F) > 0) = oo. 
But this is clearly due to trivial loops as it can be seen directly from the definition of 
that in this simple framework they cover the whole space X. 

Note however that 

li(l{F) > 0,p > 1) = > 1) - ii{T{F) = 0,p > 1) = > 1) - > 1) 

It follows that the probability no non trivial loop (i.e. a loop which is not reduced to a 
point) in intersects F equals 

eM-Mmi) > lAF) > 0})) - ' ^"^^^""^ ^« 



n.eFA.det(G)^ ' 

Recall that by Jacobi's identity, for any (n +p, n +p) invertible matrix A, 

det(A"^) det(/lij, 1 <i,j <n) = det{{A~^)k,i,n <k,l<n + p). 



In particular, det(G'^) = ^^^j^^ , so we have the 

det{G\FxF) 

Proposition 8 The probability that no non trivial loop in Ca intersects F equals 

Moreover E{e-(^--^'^)) = i^^^^^^T 

In particular, it follows that the probability no non trivial loop in visits x equals 
( a g^.^ )" which is also aconsequence of the fact that follows a negative binomial 
distribution of parameters —a and Y'k=^ 

Also, if Fi and F2 are disjoint, 

/x(/(Fi)/(F2) > 0) = >0,p>l)+ /i(/(F2) > 0,p > 1) - U F2) > 0,p > 1) 

det(G') det(G'^^^^^) 
~ det(G'^i)det(G'^2) >■ 
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Therefore the probabihty no non trivial loop in intersects Fx and equals 

exp(-.M{wo > i.nfe) > m - ( ttg^lS°? "° 

It follows that the probability no non trivial loop in La, visits two distinct points x and 
y equals ' ^I'l^^^ '"^ )" and in particular 1 — ^J^y^y if a = 1. This formula can be 
easily generalized to n disjoint sets. 

5 The Gaussian free field 
5.1 Dynkin's Isomorphism 

By a well known calculation, if X is finite, for any % e M^, 



(27r)l^l/2 J ^ ^^""^"^^ Y det(G') 

and 

(2^)|x|/2 J z z[e ^ 6 2 n^^^dz - (GJ y ^^^^^^ 

This can be easily reformulated by introducing the Gaussian field (p defined by the 
covariance £^(0^0^) = G^'^ (this reformulation cannot be dispensed with when X becomes 
infinite) 

So we have E((e-i<^''>^>) = det(/ + GMJ-^ = y/deUG^^G^ and 
E((0'^0^e~2<'?^^'X>) = (G^Y'''^ ^/det{G^G^^ Then as sums of exponentials of the form 
g-^<-,x> g^pg dense in continuous functions on the following holds: 



Theorem 9 a) The fields Ci and have the same distribution. 



b) E<^((0^(/)fF(|02)) = /E(F(£i + 7))//^'2'(d7) for any bounded functional F of a non 
negative field. 



Remcirks: 
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a) This is a version of Dynkin's isomorphism (Cf [2]). It can be extended to non 
symmetric generators (Cf [14]). 

b) An analogous result can be given when a is any positive half integer, by using real 
vector valued Gaussian field, or equivalently complex fields for integral values of a (in 
particular o; = 1). 

c) Note it implies immediately that the process 0^ is infinitely divisible. See [3] and 
its references for a converse and earher proofs of this last fact. 

5.2 Fock spaces and Wick product 

The Gaussian space V. spanned by e X} is isomorphic to the Dirichlet space 

H by the linear map mapping (ff on G^'^ which extends into an isomorphism between 
the space of square integrable functionals of the Gaussian fields and the symmetric Fock 
space obtained as the closure of the sum of all symmetric tensor powers of H (Bose second 
quantization: See [22], [18]). We have seen in theorem 9 that functionals of Li can be 
represented in this symmetric Fock space. 

In order to prepare the extension of these isomorphisms to the more difficult framework 
of continuous spaces (which can often be viewed as scaling limits of discrete spaces), 
including especially the planar Brownian motion considered in [9], we shall introduce the 
renormalized (or Wick) powers of 0. We set : (0^)" := (G'^'^)t//„(0^/VG^) where 
in the n-th Hermite polynomial (characterized by ^Hn{u) — e*"""^). It is the inverse 
image of the n-th tensor power of G^'' in the Fock space. 

Setting as before Ux — G^'^, from the relation between Hermite polynomials and 

_ ]_ 

Laguerre polynomials L„ ^ , 

E^^{x) = (-2)"n!L;^(y) 

it comes that: 

■.(rf^ := 2"n!P„-"((^)) 

More generally, if (^i, (/)2...0fc are k independent copies of the free field, we can define 
: 11^=1 ■ = 11^=1 ■ 'Pj' '■■ Then it comes that: 

1 ni+..+nfc=n j=l 



21 



Prom the generating function of the polynomials P, 



2,0- 

n , 



Therefore 



1 ni+..+nk=n " j=l 



rf''(^) = ^:(E*r^ (17) 



Note that : 0^ := Yli 0j ~ ^ These variables are orthogonal in L^. Let /"^ = — a be 
the centered occupation field. Note that an equivalent formulation of theorem 9 is that 
the fields | : (f)^ : and Ck have the same law. 

Let us now consider the relation of higher Wick powers with self intersection local 
times. 

Recall that the renormalized n-th self intersections field jCi'^ = P"''^{Ca ) = Qn'^{jC,a ) 
have been defined by orthonormalization in of the powers of the occupation time. 
Then comes the 

Proposition 10 The fields C'C and : Yli ^?)" • ^^^^ the same law. 

2 

This follows directly from (17). 
Remark 11 As a consequence, it can be shown that: 

r 

where Sk-^,k2,...k- is the set of permutations a of k = ^2 such that 

<^({Er^ + 1' ■■■ Ei"^ + ^j} ^ {Sr^ + 1' ■■■ + empty for all j. 

The identity follows from Wick's theorem when a is a half integer, then extends to all 
a since both members are polynomials in a. The condition on a indicates that no pairing 
is allowed inside the same Wick power. 
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6 Energy variation and currents 



The loop measure /i depends on the energy e which is defined by the free parameters 
C, K. It will sometimes be denoted yUg. We shall denote Zf. the determinant det(G) = 
det(MA - C)-\ Then > 0) = log(Ze) + Elog(A.). 
is called the partition function of La- 
The following result is suggested by an analogy with quantum field theory (Cf [5]). 

Proposition 12 ^ = P// 

ii) ifc,,y > 0, A ^ _r,^^^ ^^th n,y{i) = {F+fy) - ^{i) - ^{l). 

Note that the formula i) would be a direct consequence of the Dynkin isomorphism if 

we considered only sets defined by the occupation field. 

Proof. Recall that by formula (7): fj,*{p = 1,^ = x,t e dt) = e~^^*y and — 

k, 6 = X,, n e du) = I u.,y c^r n. A-^^ n.ez/pz ^.e-'^^'-du 

Moreover we have Cx,y = Cy^x = KPy and Xx = i^x + Y^y Cx,y 
The two formulas follow by elementary calculation. ■ 
Recall that //(P) = G^'^ and n{Nx,y) = G^'^C^.y. 
So we have ^l{Tx,y) = G^'^ + G^'y - 2G^'2'. 

Then, the above proposition allows to compute all moments of T and I relative to ji^ (they 
could be called Schwinger functions). The above proposition gives the infinitesimal form 
of the following formula. 

Proposition 13 Consider another energy form e' defined on the same graph. Then we 
have the following identity: 

^ _ gEAfx..log(5^)-E(Ai-Axr 

Consequently 

^^((eS^^..io^(fe)-S(^^-A^)'^ - 1)) = log(^) (18) 

Ze 

Proof. The first formula is a straightforward consequence of (7). The proof of (18) 
goes by evaluating separately the contribution of trivial loops, which equals Y^T^ogih-). 
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Indeed, 

;..((e^-^'^^-(&)-^(^-^^)'^ - 1)) = > 1) - > 1) 

The difference of tlie first two terms equals log{Ze/)+Y^ log(A^) — (log(Ze)— X] log(Aj;)). 
Tfie last term equals Ylx /o°°(^ — l)^dt which can be computed as before: 

A^e(l{p=i}(e^(^-^^)'^ - 1)) = - Elog(^) (19) 



Remark 14 (h-transforms) Note that if C'^y = h^hyCx,y and k'^ = —hLhX for some 
positive function h on E such that Lh < 0, as X' = h^X and [P']y = ^Pyh^, we have 



Remark 15 Note also that [^]^ — E(e al^' el(<A))^ is the Gaussian free field associ- 



■Z. 

ated with e. 



Integrating out the holding times, formula (18) can be written equivalently 

r" \ 7 

-J/ H AT T r r /\:r T ^^ I 1 ^ N -, / 



and therefore 



/^e( n i^f-' u^v^""'^' - 1) = log(^) (20) 

{x,y) X ^ 



Note also that Yl^xJ^^-^ = n{.,,}[g:^]"^-^+^-^- 
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Remcirk 16 These ^ determine, when e' varies with ^ < 1 and ^ — 1, the Laplace 
transform of the distribution of the traversal numbers of non oriented links N^^y + Ny^^ . 

Other variables of interest on the loop space are associated with elements of the space 
A~ of odd functions uj on oriented links : o;^'^ = —uj^'^. Let us mention a few elementary 
results. 

The operator [P*''^-']^ = exp(icc;^'^) is also self adjoint in P^(A). The associated loop 
variable writes ^^y^^'^-^^^yi^)- denote it f^u. Note it is invariant if o;^'^ is 

replaced by ui^'^ + — for some g. Set [G*^"^)]^'^ = ^ ^ ^ — By an argument similar 
to the one given above for the occupation field, we have: 

P* .j,(e*-^!'^ — 1) = exp(t(P*^'^^ — exp(t(P — Integrating in i after expanding, 
we get from the definition of fi: 

r °° 1 

/ {e'f'- - l)d„{l) = J2 Tm{P^-^f) - TriiPf)] 

k=l 

Hence 

J {e'^i^ - l)dfx{l) = log[det(-L(/ - pH)-i] 
Hence / (e'fi'' - = log[det(-L(7 - pH)-i] and 

J ju) - l)n{dl) = log(det(G'('")G'-^)) 

We can now extend the previous results (18) and (20) to obtain, setting det(G*^'^^) = Z^,^^ 

^^(g-EAf...log(g|^)-E(AL-A.)r.+»/,u. _ ^ log(^) (21) 

and ^ 

x,y ^^'f 

Let us now introduce a new 

Definition 17 We say that sets Aj of non trivial loops are equivalent when the associated 
occupation fields are equal and when the total traversal numbers X^^^a ^x,y{l) o,re equal 
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for all oriented edges {x,y). Equivalence classes will be called loop networks on the graph. 
We denote A the loop network defined by A. 

Similarly, a set L of nan trivial discrete loops defines a discrete network characterized 
by the total traversal numbers. 



Note that these expectations determine the distribution of the network defined by 
the loop ensemble Ca- We will denote B^''^'''^ the variables 

/^/ 

x,y ^^'y 

Remcirk 18 This last formula applies to the calculation of loop indices: If we have for 
example a simple random walk on an oriented planar graph, and if z' is a point of the 
dual graph X' , cuz' can be chosen such that cuz' is the winding number of the loop around 
a given point z' of the dual graph X' . Then e^'^^'e^a ■^''^^ is a spin system of interest. We 
then get for example that 

liij^ujz, 7^ 0) = J^\og{det{G^'-^-^'^G-'))du 



and hence 



leCo 

Conditional distributions of the occupation field with respect to values of the winding num- 
ber can also be obtained. 



7 Loop erasure and spanning trees. 

Recall that an oriented link g is a pair of points {g~,g^) such that Cg = Cg-^g+ ^ 0. 
Define -g = {g+,g^). 

Let ii^^y be the measure induced by C on discrete self-avoiding paths between x and 
y: fj,_l^(^x,X2, ■■■,Xn^i,y) — Cx^x'2C!x\,x3---C^xn-i,y 

Another way to defined a measure on discrete self avoiding paths from x to y is loop 
erasure (see [7] ,[19] and [8]). In this context, the loops can be trivial as they correspond 
to a single holding times, and loop erasure produces a discrete path without holding times. 

We have the following: 
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Proposition 19 The image of //^'^ by the loop erasure map 7 — > 7^^ is /i^^ defined 
on self avoiding paths by i^BE{r)) = /^^'^(^) di?^§F) = /^^'^(^) det(G|{^}x{,,}) (Here {77} 
denotes the set of points in the path r) ) 

Proof. Ur] = {xi = x,X2, ...Xn = 2/),and rjm = {x, ...Xm), 

y k=2 

where A*|^j.^ denotes the bridge measure for the Markov chain killed as it hits x and 6 the 
natural shift on discrete paths. By recurrence, this clearly equals 



as 



r^{..-.}=i.. ^ det([(7-P]|{,,}c,{,,}c) ^ det(y{''--^}-) ^ ^^{G^;;--^ 
^ det([(7-P]|{,_,}c,{,_,}c) det(y{'^'«}=) det(G{^"'}^) 

for all m < n — 1. ■ 

Also, by Feynman-Kac formula, for any self-avoiding path 77: 

/ e-<^'>^>l{,..=,}/.-'nci7) = /'^ImI /^^'^I^) = det(G,)K,}xW/^^/(r7) 
J det(Gi ) 

^ det(Gx)lWxM (-^) 
det(G'|{^}x{^}) 

Therefore, recalling that by the results of section 4.3 conditionally to ?], Ci and i2|^''^ 
are independent, we see that under /x^'^, the conditional distribution of 7 given 7^^ = 77 
is the distribution of Ci — i.e. the occupation field of the loops of Ci which intersect 

V- 

More generally, it can be shown that 

Proposition 20 The conditional distribution of the network defined by the loops of 
7, given that 'y^^ = rj, is identical to the distribution of the network defined by Ci/ 
i.e. the loops of Ci which intersect r]. 
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Proof. Recall the notation — det(G). First an elementary calculation using (7) 
shows that [i^J^ ^'^'^Is^^be^^-^ equals 



n-i,y ' 



x,X2 xi,x:j" Xn—^ 'II r , . — / i r,'-^^, 

Cx,X2Cxi,Xz---C'^^_^,y 



(Note the term e ^'^) can be replaced by n«(xr)^"^^^)- 

Moreover, by the proof of the previous proposition, applied to the Markov chain defined 
by e' perturbed by a;, we have also 



Therefore 



Moreover, by (21) and the properties of the Poisson processes. 
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It follows that the joint distribution of the traversal numbers and the occupation field 
are identical for the set of erased loops and Ci/dC^ . ■ 

Similarly one can define the image of by BE which is given by 

^%E{'n) — C'a;i,a;2---C'a;„_i,j;„'«a;„ det(G|{^}x{T)}), 

for T] = (xi, and get the same results. 

Wilson's algorithm (see [16]) iterates this construction, starting with x's in arbitrary 
order. Each step of the algorithm reproduces the first step except it stops when it hits 
the already constructed tree of self avoiding paths. It provides a construction of a random 
spanning tree. Its law is a probability measure P|y on the set STx,a of spanning trees of 
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X rooted at the cemetery point A defined by the energy e. The weight attached to each 
oriented Hnk g — {x,y) oi X x X is the conductance and the weight attached to the link 
(x, A) is we can also denote by Ca;,A- As the determinants simplify, the probability of 
a tree T is given by a simple formula: 

P5t(T) = n ^« (22) 

It is clearly independent of the ordering chosen initially. Now note that, since we get a 
probability 

reSTx,A {x,y)er x,{x,A)er 

or equivalently 

E n ^".^ n p'-=Yr\Y 

reSTx,A {x,y)er x,{x,A)er i-i-xex ^ « 

Then, it comes that, for any e' for which conductances (including k') are positive only on 
links of e. 




n n px \ n \> z, 

Ax,y)er ^y x,ixA)er I ^^^^^ 



and 



n p n ^ =§7 (23) 

Note also that in the case of a graph (i.e. when all conductances are equal to 1), 
all spanning trees have the same probability. The expression of their cardinal as the 
determinant Ze is Cayley's theorem (see for exemple [16]). 

Corollary 21 The network defined by the random set of loops Cw constructed in this 
algorithm is independent of the random spanning tree, and independent of the ordering. 
It has the same distribution as the network defined by the loops of L\ . 

This result follows easily from proposition 20. 
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8 Decompositions 



Note first that with the energy e, we can associate a rescaled Markov chain Xt in which 
holding times at any point x are exponential times of parameters A^;: Xt = x^-^ with 
Tt — inf(s, ^du = t). For the rescaled Markov chain, local times coincide with 
the time spent in a point and the duality measure is simply the counting measure. The 
Markov loops can be rescaled as well and we did it in fact already when we introduced 
pointed loops. More generally we may introduce different holding times parameters but 
it would be essentially useless as the random variables we are interested into are intrinsic, 
i.e. depend only on e. 

If D (Z X and we set F = D'^, the orthogonal decomposition of the energy e(/, /) = 
e(/) into e^{f—H^f)+e{H^f) leads to the decomposition of the Gaussian field mentioned 
above and also to a decomposition of the rescaled Markov chain into the rescaled Markov 
chain killed at the exit of D and the trace of the rescaled Markov chain on F, i.e. xj^^ — 
XgF, with S[ — inf (s, lF{xu)du — t). 

Proposition 22 The trace of the rescaled Markov chain on F is the rescaled Markov 
chain defined by the energy functional e^^^{f) = e{H^f) , for which 

a,beD 

^i^^ = — ^ Cx,aCb,x[G^]"''^ 
a,b&D 

and 

Zg = Z^dZ^{f} 

Proof. For the second assertion, note first that for any y E F, 

[H^Ty = 1.=. + M^) J2iG''r''a,y. 

beD 

Moreover, e{H^f) = e(/, H^f) and therefore 

AF = e{^>(l{.}) = e(l{.}, H%x}) = A. - E ^xAH^Tx = A.(l - pP) 

aeD 
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where pf^ = Ea,6eD = EaeD^a[^^]S is the probability that the Markov 

chain starting at x will return to x after an excursion in D. 
Then for distinct x and y in F, 

— Cx,y + ^C3;,a[if^]^ = Cx^y + ^ Cx,aCb,2/[G'^]"'''. 
a a,bED 

r nil 

Note that the graph defined on F by the non vanishing conductances Cx,y has in 
general more edges than the restiction to F of the original graph. 

For the third assertion, note also that G^^^ is the restriction of G to F as for all x,y E 
F, e^^^{G5y\F, l{a;}) — ^{GSy, [H^l[x}]) — i{x=y}- Hence the determinant decomposition 
already used in section 4.3 yields the final formula. The cases where F has one point was 
already treated in section 4.3. 

Finally, for the first assertion note the transition matrix [P^^-*^]^ can be computed 
directly and equals 

Py+ Ea,beDPaPyV^'^^"^]t = Py+ Ea,6ei) C6,2/[G'^^^^>]"'^ It Can be dcCOmpOScd 

according whether the jump to y occurs from x or from D and the number of excursions 
from X to x: 

oo 

[p^^^Ty - E( E PaiV'rbP^nPy^ + E p^i^^^rbP;) 

k=0 a,beD a,beD 
oo 

= E(E p:[GVCb,.r{Py^+ y1 p:[GVCb,y). 

k=0 a,beD a,beD 

The expansion of in geometric series yields the exactly the same result. 

Finally, remark that the holding times of x] at any point x E F are sums of a 
random number of independent holding times of Xt- This random integer counts the 
excursions from x to x performed by the chain xt during the holding time of . It 
follows a geometric distribution of parameter 1 — pi^^ . Therefore, -Arr = t-tt — r is the 

expectation of the holding times of xj^^ at x. ■ 

If X is carried by D and if we set e^^ = e + || ||l2(x) ^^'^ denote [e^j^^J^ by e^^'^^ we 
have 

= Cx,y + Yl CxMG^r, Pf = Yl PaKr''C,,x 

a,b a,bED 
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and Af'^> = A,(l-pf'^>). 

More generally, if e# is such that C# = C on F x F, and A = A* on F we have: 

a,b a,b&D 

and xP'^ = X.{l-pP'\ 

A loop in X which hits F can be decomposed into a loop l^^^ in F and its excursions 
in D which may come back to their starting point. Let /x^ denote the bridge measure 
(with mass [G^Y'^) associated with e^. 
Set 

1 oo ^ 

G^,y a,b£D n=l ^xPx a,beD 

and = ^[6, + EnM^^Px^]- 
Note that pf{l) = u^^yil) = z/f (1) = 1. 

A loop / can be decomposed into its restriction l^^^ = (^i,rj) in F (possibly a one 
point loop), a family of excursions 7^,,^j+i attached to the jumps of l^^^ and systems of 
i.i.d. excursions (7^., h < n^^) attached to the points of l^^^ . Note the set of excursions 
can be empty. 

We get a decomposition of ji into its restriction jjP to loops in D (associated to the 
process killed at the exit of D) , the loop measure jjS^^ defined on loops of F by the trace 
of the Markov chain on F, probability measures v^y on excursions in D indexed by pairs 
of points in F and measures on excursions in D indexed by points of F. Moreover, 
the integers n^- follow a Poisson distribution of parameter A^. and the conditional 
distribution of the rescaled holding times in before each excursion 7^. is the distribution 
(3n^-,T* of the increments of a uniform sample of n^^ points in [0 r^] put in increasing order. 
We denote these holding times by Ti^h and set / = A{1^^\ (7Si,{i+i)) {^iiiliii^i,h))- 

Then p — p^ is the image measure by A of 

The Poisson process £a^^ — G has intensity /x^^^ and is independent of 
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Note that is the restriction of Ca to F. 
In particular, if x is a measure carried by D, we have: 



n 



x,yi=F 



(recall that jCi^^ is the restriction of to F). Also, if we condition on the set of discrete 
loops P4^> 



where the last exponent + I is obtained by taking into account the loops which have 
a trivial trace on F (see formula (19)). 

More generally we can show in the same way the following 

Proposition 23 // C* = C on F x F, and X = on F , we denote B^''^* the multi- 
plicative functional J|[-^]^^'3'e"^»'eo'4Aa= -Aa,)_ 



Then, 



x,y ^^'y 



"^^^ x,y&F ^x,y xeF 

and 



These decomposition and conditional expectation formulas extend to include a current 
(J. Note that e^^^ will depend on uj unless it is closed (i.e. vanish on every loop) in D. In 
particular, it allows to define cu^ such that: 

The previous proposition implies the following Markov property: 
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Remcirk 24 If D — Di U D2 with Di and D2 stongly disconnected, (i.e. such that for 
any (x, y, z) & Di x D2 x F, Cx,y and Cx,zCy,z vanish), the restrictions of the network Ca 
to DiU F and D2U F are independent conditionally to the restriction of Ca to F. 

Proof. It follows from the fact that as Di and D2 are disconnected, any excursion 
measure or from F into D — Di U D2 is an excursion measure either in Di or in 
D2. m 



Branching processes with immigration An interesting example can be given after 
extending slightly the scope of the theory to countable transient symmetric Markov chains: 
We can take X = N — {0}, C„_„+i = 1 for all n > 1 and Ki — 1 and P to be the transfer 
matrix of the simple symmetric random walk killed at 0. 

Then we can apply the previous considerations to check that JC^ is a branching process 
with immigration. 

The immigration at level n comes from the loops whose infimum is n and the branching 
from the excursions of the loops existing at level n to level n + 1. Set F„ = {1, 2...n} and 
D = F" 

The immigration law (on 1R+) is a Gamma distribution r{a, G^'^). It is the law of JC^, 
and also of for all n > 1. Prom the above calculations of conditional expectations, 

we get that for any positive parameter 7, 

Prom this formula, it is clear that iO^ is a Markov process. To be more precise, note 
that for any n, m > 0, = 2{n A m) and An = 2, that G^f^ = C^'" - G^'^tG^'^" so that 
^\'si ~ T+^^^d thai for any n > 0, the restriction of the Markov chain to D„ is isomorphic 
to the original Markov chain. Then it comes that for all n, pn = 

^{F„,75„+i} ^ _i_ ^{F„,7<5„+i} ^ 2gi ^^^^ Laplacc exponent of the convolution 

semigroup Ut defining the branching mechanism equals = /(I — e~'^^)e^*(is. It is 
the semigroup of a compound Poisson process whose Levy measure is exponential. The 
conditional law of given is the convolution of the immigration law T{a, 1) with 
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Alternatively, we can consider the integer valed process Ar„(£^^">) + 1 which IS a 
Galton Watson process with immigration. In our exemple, we find the reproduction law 
7r(n) = 2~"~^for all n > (critical binary branching). 

If we consider the occupation field defined by the loops going through 1, we get a 
branching process without immigration: it is the classical relation between random walks 
local times and branching processes. 

9 The case of general Markov processes 

We now explain briefly how some of the above results will be extended to a symmetric 
Markov process on an infinite space X. The construction of the loop measure as well 
as a lot of computations can be performed quite generally, using Markov processes or 
Dirichlet space theory (Cf for example [4]). It works as soon as the bridge or excursion 
measures P^'^ can be properly defined. The semigroup should have a locally integrable 
kernel Pt{x,y). 

Let us consider more closely the occupation field /. The extension is rather straight- 
forward when points are not polar. We can start with a Dirichlet space of continuous 
functions and a measure m such that there is a mass gap. Let Pt the associated Feller 
semigroup. Then the Green function is well defined as the mutual energy of the Dirac 
measures 5^ and 5y which have finite energy. It is the covariance function of a Gaussian 
free field which will be associated to the field C\ of local times of the Poisson process 

2 

of random loops whose intensity is given by the loop measure defined by the semigroup 
Pt- This will apply to examples related to one dimensional Brownian motion or to Markov 
chains on countable spaces. 

When we consider Brownian motion on the half line, we get a continuous branching 
process with immigration, as in the discrete case. 

When points are polar, one needs to be more careful. We will consider only the case 
of the two and three dimensional Brownian motion in a bounded domain D killed at the 
boundary, i.e. associated with the classical energy with Dirichlet boundary condition. 
The Green function does not induce a trace class operator but it is still Hilbert-Schmidt 
which allows to define renormalized determinants det2 (Cf [21]). 

If A is a symmetric Hilbert Schmidt operator, det2(/ + A) is defined as J|(l + \i)e~^^ 
where Aj are the eigenvalues of A. 
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The Gaussian field (called free field) whose covariance function is the Green function 
is now a generalized field: Generalized fields are not defined pointwise but have to be 
smeared by a test function /. Still (f){f ) is often denoted J (j){x)f{x)dx. 

Wick powers :</>": of the free field can be defined as generalized field by approximation 
as soon as the 2n-th power of the Green function, |/)^" is locally integrable (Cf [22]). 
This is the case for all n for Brownian motion in dimension two, as the Green function 
has only a logarithmic singularity on the diagonal, and for n = 2 in dimension three as 
the singularity is of the order of jj^c^- More precisely, taking for example 7i^{dy) to be 
the normalized area measure on the sphere of radius e around x, 0(7rf ) is a Gaussian field 
with covariance cr^ — J G{z, z')T^^{dz)T^'^{dz'). Its Wick powers are defined with Hermite 
polynomials as we did previously: 

: 0(7r^)"' := {a'^)^ Hn C^^^ ). Then one can see that, J f(x) : 0(7r|)" : dx converges in 

for any bounded continuous function / with compact support towards a limit called 
the n-th Wick power of the free field evaluated on / and denoted : 0" : (/). Moreover, 
E(: r ■■ if) ■■ r ■■ {h)) = jG'^{x,y)f{x)h{y)dxdy. 

In these cases, we can extend the statement of theorem 9 to the renormalized occupa- 
tion field £i and the Wick square : 0^ : of the free field. 

2 

Let us explain this in more details in the Brownian motion case. Let D be an open 
subset of M'^ such that the Brownian motion killed at the boundary of D is transient and 
has a Green function. Let pt{x,y) be its transition density and G{x,y) = pt{x,y)dt 
the associated Green function. The loop measure ii was defined in [9] as 

where P^'^ denotes the (non normalized) excursion measure of duration t such that if 
0<ti< ...th < t, 

^'"'{^{ti) e dxi,...,^{th) e dxh) = pt^{x,xi)pt^-tAxi,X2) pt-t^{xh,x)dxi...dxh 

(the mass of P^'^ is pt{x, x)). Note that is a priori defined on based loops but it is easily 
seen to be shift-invariant. 

For any loop / indexed by [0 T(/)], define the measure / = Jq 6i(^s)ds: for any Borel 
set A, l[A) = Jq lA{ls)ds. As before, we have the following: 
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Lemma 25 For any non negative function f, 

/n 
G{xi, X2)f{x2)G{x2, X3)f{x3)...G{Xn, Xi)f{xi) J]^ dXi 
1 

One can define in a similar way the analogous of multiple local times, and get for their 
integrals with respect to a formula analogous to the one obtained in the discrete case. 

Let G denote the operator on L'^{D,dx) defined by G. Let / be a non negative 
continuous function with compact support in D. 

Note that is //-integrable only in dimension one as then, G is locally trace class. 

In that case, using for all x an approximation of the Dirac measure at x, local times 
can be defined in such a way that (^l, f^ = f f{x)dx. 

(l, is /i-square integrable in dimensions one, two and three, as G is Hilbcrt-Schmidt 
if D is bounded, since J J^^^G{x,yydxdy < oo, and otherwise locally Hilbert-Schmidt. 

N.B.: Considering distributions x such that / J {G{x,yyx{dx)x{dy) < oo, we could 
see that can be defined by approximation as a square integrable variable and 

l^i(tx) ) = J{G{x,yfx{dx)x{dy). 



Let 2; be a complex number such that Re{z) > 0. 

Note also that e + z (l, f \ — 1 is bounded by //, /\ and expands as an 



JV+l 



- {N+iy. 



alternating series S(- /))^ with e-^(^^) - 1 - 5(- ^T; 
Then, for \z\ small enough., it follows from the above lemma that 

/.(e-^(^^) + ^ /) - 1) = ' Tr{-{M^GM^r) 

As M^GM^ is Hilbert-Schmidt det2(/ + zM^GM^) is well defined and the second 
member writes -log(det2(/ + zM^GM^)). 
Then the identity 

^(e-(''/) + Jlf)-l)^- log(det 2(/ + zM^GM^)). 
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extends, as both sides are analytic as locally uniform limits of analytic functions, to all 
complex values with positive real part. 

The renormalized occupation field is defined as the compensated sum of all / 
in La (formally, Ca = Ca — J J^^''^ 6i^dsfj.{dl)) By a standard argument used for the 
construction of Levy processes, 

Ca: f) = lim( Y] {l{T>e} / f{ls)ds) - an{l{T>e} / f{ls)ds)) 

(we can denote lim^^o \^a,s-, //) which converges a.s. and in L^, as 



E((5](1{T>.} f f{is)ds)-aix{l{T>e} r f{ls)ds)f)^a /(1{t>.} f f{is)dsfiJi{dl) 

2 



and E(^>C„, ) = Tr{{M^GM^Y). Note that if we fix /, a can be considered as a 
time parameter and (^Ca,e, as Levy processes with discrete positive jumps approximat- 
ing a Levy process with positive jumps (^Ca,f^- The Levy exponent fi{l^T>e}{^~^''^'^ + 
^/,/^ — 1)) of ^£q,^£,/^) converges towards the Levy exponent of which is 

M(e-(^^) + (/; /)-!)). 

and, from the identity E(e"(^"'^)) = e-°'*(^"^''^^ + (''^)-i), we get the 

Theorem 26 Assume d < 3. Denoting L^, the compensated sum of all I in C^, we have 
E(e-(^«'^)) = det2(/ + M^GM^))-^ 

Moreover e"^'^"-^'-^) converges a.s. and in towards e~^^°"-^\ 
Considering distributions of finite energy x (i-©- such that J{G{x,y)^xidx)x{dy) < oo), 
we can see that (^Ca,X^ can be defined by approximation as lim;^_»oo(^>Ca, AG^x^) and 

E{(Z,xy) ^o^j {G{x,y)fx{dx)x{dy). 
Specializing to o; = |, A; being any positive integer we have: 
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Corollary 27 The renormalized occupation field Ck and the Wick square | : Yli ^1 '■ 
have the same distribution. 

If © is a conformal map from D onto Q{D), it follows from the conformal invariance of 
the Brownian trajectories that a similar property holds for the bBrownian"loop soup"(Cf 
[9]). More precisely, if c{x) = Jacohian^iQ) and, given a loop /, if T'^(/) denotes the 
reparametrized loop Ir^, with Jq" c{lu)du = s, QT'^{Ca) is the Brownian loop soup of 
intensity parameter a on Q{D). Then we have the following: 

Proposition 28 Q{cCa) is the renormalized occupation field on Q{D). 

Proof. We have to show that the compensated sum is the same if we perform it after 

or before the time change. For this it is enough to check that 

H[^iUTT>n}UT<6} / f{-is)ds-a {l{TT>v}kT<e} / f{ls)ds)ii{d-i)f) 

= a j {hrT>ri}hT<e} f {%)dsf ll{d-f) 

ikT>e}'i-TT<v f{ls)ds-a {l{T>e}UT<v f{ls)ds)n{d-f)]^) 

Jo J Jo 

l{T>£}lrr<)? / f{.ls)dsfii{d^) 
Jo 



and 



a 



converge to zero as e and rj go to zero. It follows from the fact that: 



j[l{T<e} f{ls)dsflJi{dl) 



and 



[1tt<?? 

/ f{-is)dsYi,{d-i) 







converge to 0. The second follows easily from the first if c is bounded away from zero. 
We can always consider the " loop soups" in an increasing sequence of relatively compact 
open subsets of D to reduce the general case to that situation. ■ 
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As in the discrete case (see corollary 6), we can compute product expectations. In 
dimensions one and two, for fj continuous functions with compact support in D: 

n{^aJi)-{CaJk))^ j Perl{G{xux,n)A <l.m <k)\[fj{x^)dxj (24) 

10 Renormalized powers 

In dimension one, powers of the occupation field can be viewed as integrated self inter- 
section local times. In dimension two, renormalized powers of the occupation field, also 
called renormalized self intersections local times can be defined as follows: 

Theorem 29 Assume d = 2. Let Ti^{dy) he the normalized arclength on the circle of 
radius e around x, and set = J G{y, z)n''^{dy)'n^{dz) . Then, f f{x)Q'^'^^ {(^Ca,'^£^)dx 
converges in for any bounded continuous function f with compact support towards a 
limit denoted (^C^, and 

f) (^a, h)) = / G'\x, y)f{x)h{y)dxdy. 



Proof. The idea of the proof can be understood by trying to prove that 
^ii I f{^)Qk''^'' {{^^a,''^e^)dxy) remains bounded as e decreses to zero. The idea is 
to expand this expression in terms of sums of integrals of product of Green functions 
and check that the combinatorial identities (15) imply the cancelation of the logarithmic 
divergences. 

This is done by showing (as done below in the proof of the theorem) one can modify 
slightly the products of Green functions appearing in E((5^'^^ {(^J^a, '^e^)Q'k"^ {(j^a^ ^rf^)) 
to replace them by products of the form G{x,yy{a^ya'^)^ . The cancelation of terms 
containing erf and/or then follows directly from the combinatorial indentities. 

Let us now prove the theorem. Consider first, for any xiX2--.Xn, s small enough and 
e < £i, ...En < 2e, with £j = ej if Xi — Xj, an expression of the form: 

n G(x,_^,x,)(a^,X^- j G(yi,y2)...G(y„,yiXnrfyi)-<:fc) 
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in which we define mj as sup{h,Xi^h — ^i)- 

In the integral term, we first replace progressively G{yi-i,yi) by G{xi-i,Xi) whenever 
7^ Xi, using triangle, then Schwartz inequality, to get an upper bound of the absolute 
value of the difference made by this substitution in terms of a sum A' of expressions of 
the form 

The expression obtained after these substitutions can be written 

Yl G{xi-i,Xi) J G{yi,y2)...G{ymi.^,ymiK:{dyi)...7r^:{dym,) 

and we see the integral terms could be replaced by (erf')™' if G was translation invariant. 

But as the distance between x and y tends to 0, G{x,y) is equivalent to Go{x,y) ~ 
Mog(||a; — y||) and moreover, G{x,y) = Gq^xaj) — H^" {x,dz)GQ{z,y), H^" denoting the 
Poisson kernel on the boundary of D. As our points lie in a compact inside D, it follows 
that for some constant C, for \\yi — x\\ < e, \J {G{yi,y2)7T^{dy2) — | < Ce. 
Hence, the difference A" between W and Ylixi-u^xi G{xi-i,Xi){a^^)™'^ can be bounded by 
eW, where W is an expression similar to 1^.. 

To get a good upper bound on A, using the previous observations, by repeated ap- 
phcations of Holder inequality, it is enough to show that for e small enough , C and C 
denoting various constants: 

1) J{Giy^,y2)-G{x^,X2r7^^,lidy^)7^p^{dy2) 

<C{el + {G{xi, X2f + l0g(£)2)l{||^^_^^||<^}) 

2) jG{y,,y2)'n^{dy^)n^{dy2)<C\\og{e)\' 

3) fG{y„y2)'7r^l{dy,)7r-^^{dy2) < C \\og{e)\' 

As the main contributions come from the singularities of G, they follow from the 
following simple inequalities: 

1') 

J \log{s^ + 2Re cos(^) + R^) - log(i?) f d9 

|log((£/it:)^ + 2{e/R) cos{e) + 1) I' d9 < C{{el{n>v-e}} + V(it:/£)l{R<vi}}) 
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(considering separately the cases where -| is large or small) 

2') / |log(£2(2 + 2 cos(^))) l'^ de<C |log(£) 1*^ 

3') /|log(£i cos(ei) + 62 cos(^^2) + + (£i sin(^i) + £2 sin(e2))'|' deide2 < C(|log(e)|)^ 
It can be proved by observing that for r < ei + £2, we have near the singularities 
(i.e. the values 6i{r) and 62{r) for which the expression under the log vanishes) 
to evaluate integrals bounded by C (— log(£ti))^(iM < C"(— log(£))^ for e small 
enough. 

Let us now show that for £ < £1, £2 < 2s, we have, for some integer N'^''^ 



a,T^e^/)Ql {\i-a,T^!^p) - dl,kG{x,y) ^ 



< Clog(£)^'.'=(Vi + G'(x,y)^'=l{ll,_,ll<^) (25) 

Indeed, developing the polynomials and using formula (24) we can express this expec- 
tation as a hnear combination of integrals under ]^7rf^('ia;i) j^7r^2('i|/i) of products of 

i 3 

G{xi, Uii), G{xi, Xj) and G{yj, Uji) as we did in the discrete case. If we replace each G(xi, yj) 
by G{x, y), each G{xi, Xi') by crf^ and each G(|/j, yf) by crf^, we can use the combinatorial 

cxicx + 1)...(q; + A; — 1) 

identity (15) to get the value Si^kG{x, yY^ . Then, the above results 

allow to bound the error made by this replacement. 

The bound (25) is uniform in {x,y) only away from the diagonal as G{x,y) can be 
arbitrarily large but we conclude from it that for any bounded integrable / and h, 

y(E(g^-"((Z;,<))gr""((>£^,0)) - Wx, yr '^' • • + ^ ~ ^h fimy)dxdy 

< C"v^l0g(£)^'-'= 

(as / / G{x, yy'^l^^^^_y^^^^dxdy can be bounded by Ce^, for example). 

Taking £„ = 2~^, it is then straightforward to check that / f{x)Ql'"'' {{jC-a, 7r^^))dx is 
a Cauchy sequence in L^. The theorem follows. ■ 

Specializing to a = |, A; being any positive integer as before. Wick powers of Yl^=i ^"j 
are associated with self intersection local times of the loops. More precisely, we have: 
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Proposition 30 The renormalized self intersection local times Cl and the Wick powers 



The proof is similar to the one given in [13] and also to the proof of the above theorem, 
but simpler. It is just a calculation of the L^-norm of 



which converges to zero with s. 

Final remarks: 

a) These generalized fields have two fundamental properties: 

Firstly they are local fields (or more precisely local functionals of the field in the 
sense that their values on functions supported in an open set D depend only on the 
trace of the loops on D. 

Secondly, noting we could use different regularizations to define the action of 
a conformal transformation © on these fields is given byj/ie k-th power of the 
conformal factor c — Jacobian{Q). More precisely, Q{c'^C^) is the renormahzed 
k-th power of the occupation field in Q{D). 

b) It should be possible to derive from the above remark the existence of exponential 
moments and introduce non trivial local interactions as in the constructive field 
theory derived from the free field (Cf [22]). 

c) Let us also briefly consider currents. We will restrict our attention to the one and two 
dimensional Brownian case, X being an open subset of the line or plane. Currents 
can be defined by vector fields, with compact support. 

Then, if now we denote by the complex valued free field (its real and imaginary 
parts being two independent copies of the free field), uj and J^{(l)di^(l)—(l)du,(j))dx are 
well defined square integrable variables in dimension 1 (it can be checked easily by 
Fourier scries) . The distribution of the centered occupation field of the loop process 
"twisted" by the complex exponential cxp{J2i^Ca Ii ^^"'"i^dl'^ll^)) appears to be the 
same as the distribution of the field : (jxj) : "twisted" by the complex exponential 



2 



1 

2"n! 



: (Xli 0f )" : have the same joint distribution. 




exp{J^{<Pd^<t) - (j)d^(j))dx) (Cf[14]). 

In dimension 2, logarithmic divergences occur. 
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d) There is a lot of related investigations. The extension of the properties proved here 
in the finite framework has still to be completed, though the relation with spanning 
trees should follow from the remarkable results obtained on SLE processes, especially 
[11]. Note finally that other essential relations between SLE processes, loops and 
free fields appear in [27], [20] and [1]. 
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